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1. Introduction 

In a recent work [1], the scheme of infrared regularization (developed in its original form by 
Becher and Leutwyler [2] ) was extended to the case of explicit meson resonances interacting 
with soft pions, as e.g. the first step towards a systematic inclusion of vector mesons in 
the meson sector of Chiral Perturbation Theory (ChPT). It is the aim of this article to 
provide a further extension of infrared regularization to the situation where baryons as well 
as vector mesons are present. This can be considered as a synthesis of the results in [1] 
and [2] . Such an extension is not only of interest in itself, but can be applied to a plethora 
of observables, where vector and axial-vector mesons are known to play an important role, 
such as the electroweak form factors of the nucleon. As one example we mention the 
contribution of the np loop to the strangeness form factors of the nucleon [3,4]. 

Infrared regularization (IR) is a solution to the following problem. The presence of 
mass scales which must be considered as 'heavy' compared to the masses and momenta of 
the soft pions will in general mess up the usual power counting rules of ChPT by which 
the perturbation series of the effective theory is ordered [5-7] . This observation was first 
made when baryons were incorporated in the framework of ChPT [8] (for a recent review 
on baryon ChPT we refer to [9]). The procedure of infrared regularization separates the 
(dimensionally regularized) one-loop graphs of baryon ChPT into a part which stems from 
the soft pion contribution and a part generated from loop momenta close to the 'heavy' 
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scale. The latter portion of the loop graph, called the 'regular' part, will usually not be 
in accord with the low-energy power counting, but can always be absorbed in local terms 
derived from the effective Lagrangian. It is therefore dropped from the loop graph, and 
only the first part, called the infrared singular part of the loop integral, is kept. Though 
both the vector mesons and the baryons interact as 'heavy' particles with the pions, the 
power counting is different for the two species, at least for the kinds of Feynman graphs we 
consider in this work. There, the vector mesons appear only as internal lines with small 
momenta far from their mass shell, so that the resonance propagator is counted as 0{q^) 
(where q indicates some small momentum scale or Goldstone boson mass) , while the baryon 
propagator is counted as 0{q~^), since the baryon is pushed from its mass shell only by a 
small amount due to its interaction with the soft pions and vector mesons. In the graphs 
we treat here, only one single baryon is present, with the baryon line running through 
the diagram undergoing only soft interactions. The number of the virtual vector mesons, 
however, is not fixed. The pion propagator is counted as 0(g~^), as usual, both the pion 
momentum and the pion mass being of 0{q). Appropriate powers of q are also assigned 
to vertices from the effective Lagrangian. Finally, the measure of every d-dimensional loop 
integration is booked as 0{q'^). This counting scheme applies, of course, to tree graphs, 
but also to the infrared singular, or soft, parts of the loop graphs. The regular parts of the 
loop graphs are not guaranteed to obey the power counting rules. These general remarks 
will be exemplified in the following sections. 

Before working out the case where both baryons and vector mesons appear in a Feyn- 
man diagram, we will briefly review the scheme of infrared regularization for loop integrals 
where only one heavy scale shows up. This will not only serve to give a unified presenta- 
tion of the method, but also provide some results needed for an application of the general 
scheme to the axial form factor of the nucleon. Before starting with the presentation of 
the formalism, let us mention that the loop integrals studied in this work have also been 
treated, using a different regularization scheme which is in some respect complementary to 
the one used here, in Ref. [10]. 



2. IR regularization in the pion- nucleon system 

When only pions and nucleons are treated as explicit fields of the effective theory, the 
fundamental loop integral one has to consider is 

Here, M is the pion mass (being of chiral order 0{q)) and m is the nucleon mass. Applying 
the low-energy power counting scheme outlined in the introduction, one would assign a 
chiral order of q'^~^ to this integral. We will see in a moment that only an appropriately 
extracted low-energy part of Imb will obey this power counting requirement. 

All the other pion-nucleon loop integrals are either only trivially modified by the in- 
frared regularization scheme, or they can be derived from eq. ( |2.lD (see sec. 6 of [2]). For 
example, the scalar tadpole integral containing only the pion propagator is not modified 
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at all, as there is no 'hard momentum' structure present that could lead to a nonvanishing 
regular part of this integral. Thus we have 



and a direct calculation gives 



IR _ T 

M — -'M 



r(i-f 



lir = ^—^M'^-'. (2.2) 
(47r)2 

Note the typical structure of the d-dependent power of the pion mass. For arbitrary values 
of the dimension parameter d, I^^ is in general proportional to fractional powers of M^, 
and will even diverge in the so-called chiral limit where the quark masses rriu^md go to 
zero (so that also M ^ 0) for small enough d (we mostly consider the two-flavor case here 
but the method can be trivially extended to include also strange quarks). Such terms will 
never occur in the regular parts of the loop integrals which stem from the high-momentum 
region of the integration: those parts are always expandable in the pion mass. The baryon 
tadpole integral, e.g., is 

^~ J (27r)rf /2 _^2' 

and is trivially expandable in M^, since it has no pion mass dependence at all. Therefore 
one has I^^ = 0. These remarks may serve to explain the terminology 'infrared singular' 
vs. 'regular'. 

Returning to eq. (p.l|), we must extract the part of this integral that is proportional 
to some d-dependent power of M, like in eq. (p.2|). One can think of this extraction 
prescription as an operational definition of infrared regularization. The method is explained 
in full generality in [2]. Here, we concentrate on the case with on-shell momentum p, i.e. 
p2 = ,^2_ shows all the features we need for the demonstration, and also yields the 

result we will use in our application of the scheme in sec. |^. We introduce a Feynman 
parameter integration in the usual way: 



" ' (2^ io [((p - 0^ - + {l^ - Af2)(l - z)]2 ■ ^^-^^ 



Performing the standard steps, we find for p2 = 77^,2. 

= -„/-*n^ /' (2.4) 

(47r)2 JO [z^-az + aY 2 

where we have defined a = M'^/rnP (note the difference to ref. [2], where this letter is 
reserved for obl = M/m). Fractional powers of the small variable a will be produced near 
z = 0: there, the integrand is approximately (a) 2"^. For small enough d, there would be 
an infrared singularity for M — > located in parameter space at z = 0. It can already be 
seen from eq. ( p. 3D that the parameter region near z = is associated with the low-energy 
portion of the integral: In this region, only the 'soft' pion propagator is weighted in the 
loop integration, while the hard momentum structure of the nucleon propagator dominates 
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near z = 1. The extraction of the part of the integral proportional to d-dependent powers 
of a now proceeds as follows: the parameter integration is split into two parts like 

= - ■ (2-5) 

/O JO Jl 

We will first show that the second integral on the r.h.s. is regular, i.e. expandable, in the 
variable a. This is easy to see, because for z > 1, the integrand can be expanded like 

Interchanging integration and summation (which is a valid operation at least for some 
range of d), one gets for the regular part 

U . „/-£Pzi) r ^ -'-'^^ ^ Oia, (2.7) 

(47r) 2 Jl [z^ -az + af-2 (3 - d) (Att) 2 

At this point we should make the remark that the extension of the parameter integration 
to infinity will lead to divergences as d increases. The infrared singular or regular parts are 
then defined as follows: The parameter integrals are computed for the range of d where 
they are well-defined, and the result will be continued analytically to arbitrary values of 
d. This amounts to the suppression of power divergences of the parameter integrals, which 
have nothing to do with the infrared singularity at 2; = 0, and will be cancelled anyway on 
the r.h.s. of eq. (^^). 

Next we must show that the first term on the r.h.s. of eq. ( ^.5| ) is proportional to a 
d-dependent power of a. To see this, we substitute z = ^/ay to get 

dz /■°° dy 



lo {z^-az^af 2 Jo (y^ _ + 1)2 % 

The remaining integral on the r.h.s. will not produce d-dependent powers of a, since 
the integrand can be expanded in ^/a similar to eq. ( [2.61 ). Thus we have found that the 
parameter integral from zero to infinity is proportional to a d-dependent power of a. The 
infrared singular part of the loop integral therefore equals 



^MB — — — ~r 



(47r)2 Jo {z'^-az + af' 



2 — k\ 

^ k=0 



As d ^ 4, this has the well-known leading term 



where the dots indicate terms of higher order in M/m. This is clearly in accord with the 
low-energy power counting. In contrast to that, the first term of the expansion of the reg- 
ular part R obviously violates this counting when d ^ A. The infrared regularization now 
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prescribes to drop R from the loop contribution and substitute for Imb- Moreover, 
the poles of /^^^ in d — 4 are also absorbed in a renormalization of the masses and coupling 
constants of the effective Lagrangian. Again, for a more general and comprehensive treat- 
ment of the IR scheme in the pion-nucleon system, the reader should consult the original 
article of Becher and Leutwyler [2] . 



3. IR regularization for vector mesons and pions 

In this section, we consider another case of infrared regularization, first examined in [1]. 
The internal baryon lines from the preceding section are now replaced by vector meson lines, 
however, the vector mesons do not show up as external particles in the graphs we consider 
here (an example for the treatment of such a graph can be found in sec. 11 of [1]). The 
only external particles here are pions (or, in some cases, soft photons etc.), so that there 
are only small external momenta of order 0{q) flowing into the loop. The fundamental 
scalar loop integral is in this case 

{2^Y{[q-iY-Miw-M^y ^^-^^ 

where My is the mass of the heavy meson resonance and q is some small external momentum 
(small with respect to the resonance mass My). This is, in principle, the same function as 
in eq. (|2.1| ) , so the reader might ask why we devote this section to the examination of this 
case. The point is that the extraction of the infrared singular part of Imv must proceed 
along different lines here. As shown in [1], a splitting like that of eq. ( |2.5| ) does not amount 
to a separation into infrared singular and regular parts as in the preceding section. This 
can be traced back to the fact that the extension of the parameter integrals to infinity leads 
to a singular behaviour of the loop function near = 0. Of course, also ij^^ such 
a singularity as the external momentum squared goes to zero (see sec. 5.4 of [2]), but the 
point = lies far outside the low-energy region in that case. Here, however, the point 
= lies at the center of the low-energy region, so the infrared singular and regular parts 
can not be expressed as parameter integrals from zero or one to infinity, respectively. For 
example, a 'regular' part defined as being proportional to the parameter integral from one 
to infinity would not be expandable around q^ = 0, as it should be in order to be able to 
absorb the corresponding terms in a renormalization of the local operators in the effective 
Lagrangian. 

To circumvent this difficulty, we will take up the following simple idea from [1]. We 
observe that Imv is analytic at q^ = 0, the only singularity being the threshold branch 
point at = (My + M)^, which is far outside the range of small g'^-values. Expanding 
Imv in (f' ■, the analyticity properties in that variable are obvious. Each coefficient in this 
expansion can be split into an infrared singular and a regular part almost like in eq. (p.5|). 
Extracting the infrared singular part proportional to d-dependent powers of the pion mass 
of each coefficient, and resumming the series, one arrives at a well-defined expression for 
the infrared singular part of the loop integral Imv that is (by construction) expandable in 
the small variable q^ /My. 
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To start with the analysis, let us first consider the special case where the external 
momentum vanishes: q = 0. Then we have 



iMviO) 



dH 



dH i f dH 



In the last step of this equation, we have already achieved the splitting into an infrared 
singular and a regular part, which is quite trivial here. The first term is proportional to a 
d-dependent power of M (compare eq. (^.2|)), while the second part is clearly expandable 
in (because My S> M^). Moreover, the pion mass expansion of the first term starts 
with ikf^"^, in agreement with the low-energy power counting for Imv'- In contrast to the 
baryon propagator in sec. ^, the resonance propagator is counted as 0{q^) since, in the 
low-energy region, the off-shell momentum of the internal resonance line is far below its 
mass shell. The counting for the pion propagator and the loop measure is the same as 
before, and one is lead to the power counting result q'^~'^ for Imv- Therefore we can write 

^w(O) = ^2^M2 / ^J2^- (3-2) 

It should be clear that there will only be corrections of 0{q^) to this result when we compute 
the regularized integral for nonvanishing q. We introduce the following variables, 

a = — TT- , 3 = — ^, (3.3) 

which we assume to be small in the sense that a,f3 <C 1. Just like in sec. we use the 
Feynman parameter trick to write 

(47r)2 h (/3z2 + 2(i_a-/3)+Q)2-- 



2 



We rewrite this expression as follows: 



,_4r(2-f) dz / , \"2 



(47r)i Jo (z(l - a - /3) -Fa)2-i V z{l - a - l3) + a 
This can be expanded according to 



(^2)^_<:!g, r(f-i) n2 -i)dz f pz' 



(47r)i ^ fc!r(f -l-k)Jo {z{l -a-p)+ af~l \z{l-a-P) + a 



In the next step, we extract the infrared singular part of each term in the sum. In each 
parameter integral, we substitute z = ay and extend the integration range to infinity: 



4= / dz / dy ^ J. (3.5) 

Jo (2(i_a-/3)+Q)2+^-2 Jo (i + y(i_a_/3))2+fc-2 
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Divergences of the parameter integral due to the extension of the upper hmit to infinity 
are again handled by analytic continuation in d. It can be expressed in terms of Gamma 
functions: 

y2fc r(2A; + l)r(l - ^ - A:) 

dy- - 







(1 + y(l _ a - /3))2+^-t (1 - a - Pr'+'Ti2 - I + k) 



This is clearly expandable in the small variables d and /3, so that the r.h.s. of eq. ( |3.5| ) 
in fact has the proper form of an infrared singular contribution, being proportional to d- 
dependent powers of d. Moreover, it is not difficult to see that the parameter integrals 
from 1 to infinity are completely regular in the small variables. So, putting pieces together, 
and using the following identity for Gamma functions: 

we can sum the series of infrared singular terms and write 

riR(,2._ K~' ^ r(2fc + i)r(i -i-k) 



(3.7) 



Note that the leading term in this result is of chiral order 0{q'^~'^), as predicted by the 
power counting scheme. All the terms we have separated off from Imv ^^^6 regular in both 
small parameters, and only the infrared singular terms remain. For /? = 0, the result for 
-^MV'(O)' which was already established in eq. ( |3.2D , is reproduced. 

It is perhaps worth noting that the result of eq. ( |3.7] ) can be obtained in a different 
way, which goes back to Ellis and Tang [11, 12]. Though they use their method in the 
pion-nucleon sector, a variant of it is also applicable here. Their prescription to obtain the 
soft momentum contribution of a loop integral is the following: Expand the propagators of 
the heavy particles as if the loop momentum were small, and then interchange summation 
and loop integration. It is claimed that this prescription eliminates the 'hard momentum' 
contributions present in the full loop graph. If this is true, and the concept of hard vs. soft 
momentum effects is a well-defined one, the result of the procedure should reproduce the 
infrared singular part of the loop integral in question. Let us see how this works out for 
our example. Following the prescription just described step by step, we make the following 
set of transformations: 



{2ttY /2 _ M2 ^ (g2 + /2 _ m2 )'=+! 

d'h i ^ {2q-l)^ 



k=0 



{2TrY ^2 _ M2 ^ (g2 + M2 - m2 )fe+i 
d'^l i{2q ■ /)'= 



{2Try (/2 - M2)(g2 + M2 - M2 )fc+l ' 
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While the first step is just the expansion of the vector meson propagator pole imposed by 
the prescription, the second step deserves a comment: There, we have used the same trick 
of partial fractions to split off some hard momentum contributions as in the treatment 
of /mv(0)) but now this was performed k + 1 times. In the last step, summation and 
integration were interchanged. Using the formula 

p^ra = '-^ ) r(i) (4;)j 

and the fact that the loop integrals in the series vanish if k is odd, the result of the 
transformation is 



, 2. M'^-' ^ {-Aq'M'r r(n + l)r(l - f - n) 



Extracting a factor of My and using the identity 



r(n + i) r(i) ' ^'-''^ 



we see by comparing to eq. (3/7) that Iffy = ij^y- 



In [1], the result for ij^y was given in a different form (see eqs. (7.10) and (8.4) of that 
reference). Of course, it is equivalent to the result derived above: it is shown in app. ^ 
that the two different forms just amount to a reordering of the corresponding expansions. 
We will find that the form of eq. (3.7) is most practical for our purposes, in particular, the 



chiral expansion can almost immediately be read off from that formula. 

In closing this section, we note that the spin or the parity of the resonance obviously do 
not play a major role in the above considerations. Though we will concentrate on the case 
of vector mesons, most of what we have said would also apply for other meson resonances, 
like e.g. scalar or axial-vector mesons. 



4. Pion-nucleon system with explicit meson resonances 

Now that we have collected the results for the fundamental one-loop integrals in the pion- 
nucleon and the vector meson-pion sector, we are prepared to extend the framework of 
infrared regularization once more, and apply it to Feynman graphs where nucleons, pions as 
well as vector mesons take part in the same loop. The simplest example where this situation 
occurs is the triangle graph consisting of one pion, one nucleon and one vector meson line. 
Due to baryon number conservation, the nucleon line must run through the complete 
diagram. We shall assume that only a small momentum k (small in the usual sense) is 
transferred at the vector meson-pion vertex. Such a graph will typically contribute to some 
nucleon form factor in the region of small momentum transfer. With this application in 
mind, and with the excuse that it will simplify the presentation a bit, we will further specify 
to on-shell nucleons. Let p and p be the four-momenta of the incoming and the outgoing 
nucleon. Then we have 

p'^ =111^ =p^ = {p + kf ^ k'^ = 2p-k = -2p ■ k. (4.1) 
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Using these kinematic relations, we can rewrite the fundamental scalar loop integral 

with the help of the usual Feynman parameter trick, as 
df^l f^-y 2idxdy 



r dH r r 

J j2^Jo Jo 



{y{{p - I? - m?) + x{{k + 0^ - M2 ) + (1 - X - y){P - M2))3 ' 
Doing the loop integration in the usual manner, we get 

r(3-f) f^-y dxdy 



(4.3) 



lMBv{k ) =m ^ 
(47r; 

where the definitions 



(47r)2 Jo Jo (y2 _ ay + a + /3(x2 + xy) + x(7 - a - 2 

(4.4) 



a = ^<Cl, /? = — <1, ^ = _^^oi. 

were used. In particular, we have assumed here that the nucleon and the meson resonance 
are roughly of the same order of magnitude (in the real world, we have 7 ~ 2/3 for the rho 
resonance, which is good enough for our purposes). 

We should remark here that there is, of course, a second graph with the same topology, 
where the pion and the resonance line are interchanged (see fig. |2|). The expression for the 
corresponding scalar loop integral is 

iMBvik^) = j (2vr)d {{p - 1)2 _ m^)^(l _ k)^ _ mI){P - AP) ' ^^'^^ 
However, the reader can convince himself that, due to the on-shell kinematics specified in 



eq. (4.1), this will give exactly the same expression as in eq. ( ^.41 ). Thus, we can focus on 
the integral Imbv- 

In analogy to sec. ^ it will be instructive to begin with the special case where A: = 0. 
Since we must split off the terms where only propagators of heavy particles occur, we can 
obviously apply the same partial fraction method that led to eq. ( |3.2| ). The remaining 
pion-nucleon integral can be dealt with as in sec. S. This gives 



(0) = • (4.6) 



'MBV\-) - M2-MI 



Having the standard method described of sec. ^ in mind, we note that this equals 

I^^ (0) - m'^-'^^'' - i'^ l'^ ^^^^ 
J^MBVV^) — ^ 



(47r)2 JO Jo {y"^ — ay + a + x{j — a))" 
r(2-f) dy 



(7 -a) (47r)2 Jo {y"^ -ay + af~2 

(compare the last line with the l.h.s of eq. ( |2.8D ). Once again, the possible divergence for 
large d at x — > 00 was regularized by analytic continuation from small d as before. It is 
reassuring to see that eq. (|4.60 is reproduced in this way. 



-9- 




The extension of both parameter inte- 
grations to infinity is the natural generahza- 
tion of the prescription used in sec. ^ This 
can be seen as follows. In eq. (4.3), we chose 
our Feynman parameters such that the pion 
propagator contributes with the weight one 
at X = y = 0, while the vector meson and 
the nucleon propagator have their maximum 
weight at {x,y) = (1,0) and (0,1), respec- 
tively. Consequently, the infrared singular- 
ity is located in parameter space at the point 
{x,y) = (0,0). Indeed, a look at eq. (U) 
confirms that the integrand in that expres- 
sion is approximately a2~^ in the region 
where x ~ y ~ 0, thus producing infrared 
singular terms there. Far from this region, 
the integrand is expandable in a, as we will 
show later. To be more specific, imagine the 
positive quarter of the parameter plane, i.e. 
Q+ = {{x,y) : X > 0, y > 0}, split in two 



Figure 1: Illustration of the integration ranges 
for the infrared singular and the regular part of 
Imbv- The shaded triangle is the region inte- 
grated over in eq. (4.4). The three rays indi- 
cate the extension of the connections from the 
soft point at (0, 0) to the hard line from (0, 1) to 
(1,0). 



parts, the first being the triangle one integrates over in the full loop integral (see eq. ( |4.4| )) 
and the second part its complement in named i?+ (see fig. || for an illustration). In 
the latter parameter region, no infrared singularities are located, hence, the integration 
over this region should yield a proper regular part (there is a qualification here, see below). 
Geometrically, one might imagine that all lines from the 'soft point' (0, 0) (the location 
of the infrared singularity in parameter space) to the 'hard line' from (0,1) to (1,0) are 
extended to infinity to get the infrared singular part of the loop integral. In a more general 
setting, there may be soft or hard points, lines, surfaces etc., depending on the number of 
soft and hard pole structures in the integral under consideration. Extending all connecting 
lines from the soft point, line etc. to the hard point (line. . .) to infinity, one always achieves 
a splitting in the original parameter region and a region where no infrared singularities in 
parameter space are present, analogous to It is possible to show that this geometric 
picture leads to the same results as the prescription given in sec. 6 of [2]. This picture 
might serve as a guide when trying to split arbitrary one-loop graphs into infrared singular 
and regular parts, however, one should always convince oneself that both those parts have 
the correct properties, and that their sum equals the original integral. One has to be a 
bit careful here, as the following qualification shows. As explained in sec. ^, extending the 
integration range to infinity might lead to unphysical singularities in variables in which the 
original integral is analytic (at least in the low-energy region). Those singularities will be 
harmless if they are located far from the low-energy region, as e.g. the singularity at s = 
of the infrared singular parts in the pion-nucleon sector, but in other cases, they can be 
disturbing. Following the method of [1], we avoided this problem in sec. |^by expanding the 
original loop integral in the small variable (5 beforehand. Only then could the integration 
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range be extended to infinity, in each coefficient of tlie expansion. We must expect that a 
similar phenomenon will occur in the present case. 

To exclude this from the start, we expand Imbv in analogy to eq. 



^^^^ (47r)i j^r(f -2-i) io io j\{y2 _ay + a + x{-i -a- 0)f-i+^ 

(47r)i j'^r(f -2- j) io io j!(y2 _c,y + a + x(7-a-/3))3-f+j' 

There is still some /3-dependence in the denominator, but in that combination, it will turn 
out to be harmless. We extend the parameter integrations to (5+ and define 

riR (u2^ _ ^ r(f - 2) r(3 - i)rfxdy(Elo 



f 47r^ 2 tr; i 



(47r)2 ^t^r(f -2-j) Jo ^ j!(y2_ay + a + x(7-Q-/?))3-2+i 

Since the denominator is always positive for 0<Q!^1,|/3|<C1, the x-integration can 
readily be done using the formula 



x" _ a"+^--Pr(n + l)r(D - (n + 1)) 
{a + bx)D ~ ~lf^^ f{D) ' ^ ' 

together with eq. (^). This gives 

riR 2^ _ + i + i)r(2 - f - /) i-py 

^'^^^^"(4vr)i,t^t^ r(, - / + l)ra + 1) (7-a-/3)^+'+i ^ 



In the next step, we can write down the chiral expansion of the y-integral, using the same 
method as in sec. 0. The generalized formula is 



dy y = V ^ r(^)r( , ) 



Obviously, the chiral expansion of this integral can straightforwardly be read off from the 
series on the r.h.s. Inserting this result, we get 

^^^^^^ ^" (4^)f 2r(i-/ + i)r(fe + i)r(/ + i) " 

(4.9) 

The only expressions we have not expanded in the small variables are the factors of (7 — 
a — /?) in the denominator, but this can of course be done: the corresponding geometric 
series is absolutely convergent due to the assumption that 7 ~ 0(1). 

To complete the proof that eq. (|4.9| ) is the correct infrared singular part of 1m bv > we 
have to show that all the terms we dropped in the extraction procedure described above 
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are regular in a. Those terms are proportional to parameter integrals over the region 
of the general type 



dxdy{x{x + y)y 



R+ {y'^ - ay + a + x{j - a - (3))^ 2+i 

dxdz{xzy 



lz=i Jx=o {{z - x)2 - a{z - 1) + x(7 - 
Here we have traded the variable y for z = x + y. One finds that the function 

z — I 

has the property 

f{z,x) > min{4, (7-/?)} 

in provided that the parameters f3 and 7 are in their typical low-energy ranges. This 
is already sufficient to ensure that the integrand of Rj can safely be expanded in (z — l)a, 
and that integration and summation of the corresponding series can be interchanged. This 
proves the regularity of the integrals Rj in a. 

It is also possible to show that the result for ll^^y can be obtained in a way that is 
closely analogous to the prescription of Ellis and Tang (see the end of sec. ^). The proof 
can be found in app. 



5. Application: Axial form factor of the nucleon 

A typical example for an application where the triangle graph treated in sec. ^ shows up is 
a contribution to some form factor of the nucleon at low momentum transfer. As a specific 
example, we consider the nucleon form factor of the isovector axial-vector current in a 
theory with an explicit rho resonance field. A representation for this form factor, using the 
infrared regularization scheme in the pion-nucleon sector, has been given by Schweizer [13]. 
In that framework, the contributions due to the various baryon or meson resonances are 
contained in the low-energy coefficients (LECs) of the effective pion-nucleon Lagrangian, or, 
more correctly: The contributions from tree-level resonance exchange can be described as 
an infinite sum of contact terms derived from the pion-nucleon Lagrangian. The inclusion 
of explicit resonance fields therefore amounts to a resummation of higher order terms, which 
is often advantageous (as discussed in detail in ref. [14]). Moreover, a theory with explicit 
resonance fields can serve to achieve an understanding of the numerical values of the LECs, 
relying on the assumption that the lowest-lying resonances give the dominant contributions 
to those coefficients. This is usually called the principle of resonance saturation, and has 
been very successfull in ChPT, see e.g. [15,16]. Assuming this principle to be valid, the 
pion-nucleon LECs can be expressed through the masses of the resonances and the couplings 
of the resonances to the nucleons and pions. Such relations are most useful, of course, if the 
resonance masses and couplings are sufficiently well known. In this section, we will compute 
the contribution to the axial form factor of the nucleon that is described by the two triangle 
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A(k)> A(k) 




V(l-k) // \ P(l) P(l) / \, V(l+k) 



N(p) p-1 N(p) N(p) p-1 N(p) 

Figure 2: Triangle graplis contributing to tlie axial form factor of the nucleon. The dashed line 
represents the pion, while the double line stands for the vector meson. The letter A indicates the 
external axial source, N denotes the nucleon under consideration. 

graphs of fig. (|2|), with one nucleon, one pion and one resonance line, in the framework 
of the extended infrared regularization scheme developed in the preceding sections. First, 
we must set up the necessary formalism and collect the various terms from the effective 
Lagrangian we need for the computation. 

By Lorentz invariance, the matrix element of the axial current 

A'^{x) = q{xhf,-f5—q{x) 
between one-nucleon states can be parametrized as 

{N'{p)\^^{x)\N{p)) = u'ip) (gaWtm + Gp{t)^ + GT{t)^yj-n{p)e^'^\ (5.1) 

In the above expressions, q is the quark field spinor, q^ = {u, d), N' is the outgoing nucleon 
with momentum p, N labels the incoming nucleon with momentum p, and k = p—p, t = k"^. 
The symbols r* denote the usual Pauli matrices. Finally, u' and u are the Dirac spinors 
associated with the outgoing and incoming nucleon, respectively. Assuming perfect isospin 
symmetry and charge conjugation invariance, as we will do here, leads to Gt = 0. The 
relation of Ga and Gp to the quantities Fi^2 used in [13] is Fi{t) = GA{t),2mF2{t) = 
—Gp{t). For later reference, we give the representation of Ga up to order that can be 
found in [13]: 

2 , , . 9\M^ 



GA{t) = QA + ^IgM' + d22t 
9aM^ 



167r2F2 



Here F and qa denote the pion decay constant and the nucleon axial charge in the chiral 
limit, respectively, while the coefficients Ci,di are LECs showing up in the pion-nucleon 
effective Lagrangian at order two and three, respectively. For a precise definition of the 
underlying Lagrangian see ref. [17]. 



-13- 



We now turn back to the calculation of the axial form factor in the presence of vector 
mesons. We write down the relevant terms in the effective Lagrangian and give the neces- 
sary rules for the vertices and propagators required for the calculation. First, the lowest 
order chiral Lagrangian for the pion-nucleon interaction reads 



Here, ■0 is the nucleon spinor, ijp- = (p, n), and the matrix collects the pion fields vr" via 



£5^) = ^{ilp - m)^ + ^#75V'. (5.3) 



exp 



IT ■ TT 



2F 

In the last line, we have introduced external isovector vector and axial-vector sources, v 



and a^, 



fj. 2 ^ 2 



The covariant derivative in eq. (|5.3| ) is defined as 



^From eq. (^.3|), one derives the NNir vertex rule 

9A 
2F^ 



for an outgoing pion of momentum q. 

Now we turn to the effective Lagrangians involving the vector meson fields. We choose 
a representation in terms of an antisymmetric tensor field W^u [1,15,18]. The corresponding 
free Lagrangian is 

where 

The brackets (...) denote the trace in isospin space. From C^, one derives the tensor field 
propagator in momentum space. 



At lowest chiral order, the interaction of the rho meson with the pions is given by [15] 
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In the first term we have used the definitions 



F^U ~ dfj,lu dyl^ ^[^/J) 

F^jj — d^Ty dijT^ ^['"/^) ^v]' 

The external sources r^j_,l^ are counted as 0(g), so that is of chiral order 0{q^). Also, 
Ufj^ is of 0{q). Therefore, leads to vertices of chiral order 0{q^). Using the method of 
external sources to derive Greens functions from the generating functional, we must extract 
the amplitudes linear in the source to compute the matrix element of eq. ( |5.1| ). For the 
triangle graphs considered here, we need the vertex that connects the vector meson and 
the pion with the external axial source. From eq. ( |5.4| ), we find the corresponding vertex 
rule 

laci Gv , _ _ _ _ ^ Fy ^ 



Here i,a,c are the isospin indices associated with the axial source a-r, the pion and the 
vector meson field W^y, respectively, q is the four-momentum of the outgoing pion. Since 
k and q are counted as small momenta, the chiral order of this vertex rule is in accord with 
the power counting for the interaction Lagrangian. 

The restrictions of chiral symmetry are not that strong for the interaction of the vector 
mesons with the nucleons: Here, there are terms of chiral order 0{q^). The leading terms 
of the interaction Lagrangian have been given in ref. [19]. For the SU{2) case we consider 
here, the relevant terms are 

Cnw = Rvi"T^.uW''^^ + 5yVS7^D,W^^> 

+TviJ^^,DxW''''D^D,^ + Uvi^ax^W^''' D^D^^. (5.5) 

In the notation of [19], we have Ry = Rd + Rf,Sv = Sd + Sf, etc. The definition of 
a^i, is standard, a^u = i[y^,^^]/2. It turns out that only the piece proportional to Gy 
from eq. (|5.4| ) and the piece proportional to Ry from eq. ( |5.5| ) contribute at lowest order 
to the diagrams computed here, the chiral expansion of which starts at 0{q^) (given that 
a scheme like infrared regularization is used that preserves the power counting rules). To 
keep the presentation short, we show only the contribution from those terms and therefore 
neglect some higher order contributions. However, this will be sufficient to compare our 
results to the representation up to O(g^) given by Schweizer [13]. 
The evaluation of the first graph gives (see fig. (|2|)) 

and the second one gives 

dH Ry^_ JT^7'"'i^ + k)/gA,_A i (Gy^^^ ,^ ,,,, 



T^J^P-{l + k) fgA , a\ i (Gvn 
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We have left out the Dirac spinors u, u here. Since we consider on-shell nucleons, we 
can use the Dirac equation to simplify the numerators of the integrals, u{p)p = u{p)m, 
]j>u{p) = mu{p). We will make some remarks on the computation of Ii (the computation 
of I2 can be done analogously). In a first step, we reduce the full loop integral to a linear 
combination of scalar loop integrals, which have been treated in detail in the preceding 
sections. Scalar loop integrals without a pion propagator denominator are dropped using 
infrared regularization, since they are pure regular parts. Therefore we can replace 
everywhere in the numerator. 



With the abbreviation 



^ 2V2gAGvRv 



we get 

, d'^l 2imt{{M^ -k-l){lP -kP)apr+'^{kr-lr){t^-U))l^ 

h = 91'^ 
+ 



{2^Y Hp - If - m^M - - M2 )(/2 _ M2) 

dl^l i{{M^ - k ■ 1){IP - kP)apr + - lr){t^ - 



{2^Y {(I - kf - M^)[P - M"^) 



{2ttY {{p - ly - m^){P - AP) ^ 

For completeness, we shall give the relevant loop integrals with tensor structures in app. 
Using the coefficient functions defined there, the result for the first integral Ii can be 
written as 

h = -igiT^^rl'C^ + A;,/f ) + p./{P))75, (5.6) 
where the coefficients read 

/{^^ = 2m{Ml - k^)Ci + mk\4jsv + Imbv){M^ + - 
-mM^{M^ + M^- k^)lMBV + mM^lMB, 

/f ) = 2m^(M^ - k^)i3C2 + C3 + 4C4) + Am^M^I^jBV + Imbv) 

-2m\Ml - k^mt.BV + Imbv) + '^m\lt,BV " Imbv){M^ + - k") 
-Am^M^lMBV + 4iV + - + - M^)Imv 

+{M^ + M^-k^){lMV-lii^v)^ 
= 2m\M^ - fc2)(3C2 -C3- 2C4) + 4m^M^{I^By - ig^y) 

-2m^{li,BV - Imbv){M^ + M^- k') + 2t% - 2m2/«g 
+(m2 + mI- k^){I^j^^y - Imv) - hi. 

What concerns the evaluation of I2, we note that it is given by 

h = -i9ir*(7r4'^ + kr4'^ +Pr/?^)75, (5.7) 

with 

r(7) _ r(7) Ak) _ Ak) Ap) _ _ Ap) 
-'2 — -'1 ' -'2 — -'1 ' -'2 — -'l • 
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The sum of both graphs therefore gives 



h+2 = h+h = -igir\^^I^{i^ + krl'Cl^)^^ 



with 



(7) 
1+2 



r(7) 
'l ' 



r(fc) 
'1+2 



In the sum of the two graphs, the contribution proportional to {p + p)r cancels, as was 
to be expected on general grounds (see the remarks following eq. (^.iP). We are now in a 
position to display the decomposition of the graphs as a linear combination of the scalar 
loop integrals worked out in the preceding sections: 



7-(7) _ Si) tIR _l „{7) r/i?, , „(7) rlR 
-'1+2 — '^MBV^MBV "T (^MB^MB "r Cjvfy-'Afy) 



'1+2 



ik) jiti ir'^'^l Tin I Jin I r.^ 

^MBV^MBV "T Cmb-^MB + C^y-'MV + Cj^^ -'M 



tIR 



Sk) tIR 



(k) tIR 



Sk): 



(5.8) 
(5.9) 



The expressions for the coefficients c^'''-* read 



J7) 
'-MBV 



4m 



(d- 2)A;2(A;2 -4m2) 



m'^M^ + {k^{{d - 5)m^ + M^) - 2m^M^)M^ 



2 )\/^2^ 



((M^ + m'\2d - 7))k'' - 2{d - 2)rn^M'^k^ - m^M'')M'^ 



,2 



+ {d-^)k^m^{k'^ - M'^) 



r2\2 



.(7) 
-MB 



.(7) 



2m 



(d- 2)fc2(A;2 _4^2) 
2m 



(M^ + (d - 3)A;2)((2m2 - M^)^^ + 2m2(M2 - M^)) 



(d-2)(A;2 -4m2) 
and the coefficients c'-'^-* are given by 



{Ml + {d- 3)k^)ie - M2 - M^; 



(fe) 



2m2 



(d-2)A;4(A;2-4m2) 



{{d - 2)k'^ - 4:{d - l)m'^)Ml 



- 2{{d - 2)k^ + {dM^ - 2{d + I)m2)fc2 _ 4(d _ l)m'^M^)M^ 

+ ((d - 2)k^ + 4((i - 5)m2fc4 - 2{d - 4)M^k^ + {{d - 2)k^ - 4(d - l)m'^)M^)Ml 



A{d-2,)k^rr?{k'^ - M'^) 



r2\2 



1 



=MB (d-2)fc4(A;2-4m2) 

+ 2n?{[d - 2)M^ + {d- A)M'^MI + An?{{d - 'i)M^ + 2M^))A;2 

+ 8((i - l)m^(M2 - M^)M^ , 



(4(d - 3)m2(2m2 - M'^) - {d - 2){2m^ + M'^)Ml)k^ 



Sk) 



-M 



4m2 



(d-2)A;2(A;2 -4m2) 



(M^ + (d - 3)A;2)(A;2 - Af 2 _ ^2 ^ 



)t2 ■ 
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Here we used the abbreviations A;^ = (A;^)^ and = (A;^)^. In view of the denominators of 
the coefficients c^'^'^\ which contain powers of k"^, it is advantageous to expand the scalar 
loop integrals in the small variable f3 = 7/? first. From eqs. (|3 .71,14 .91), we find 



and 

tIR 

tIR _ ^MB 
^MBV 



_ Im , ^f{d'j-{d-A)a)lM 



m?{a — 7) 

J {{d - 1)(7 - a){a - 2) - 2(4a - a^))m^li^j, + {{d - 3)a - {d - 1)7)/m \ 
''^V 2(d-l)m4(7-a)3 J 

+ 0(/3'). 

Inserting the /^-expansions of the scalar loop integrals in the expressions for /];_[.'2 from 



eqs. ( 5^ , 5^ ), one observes that the poles in the variable /? ~ /c^ cancel. In the final step, 
we must insert the expansion of the scalar loop integral in the second small variable 
M^, which can directly be read off from eq. ( |2.8D . The expression for Im is given in 
eq. ( |2.2| ). Doing this, taking the limit d — > 4 and comparing to the decomposition of the 
matrix element in eq. (^.1]), one finds the following 0(g^)-contribution of Ii and I2 to the 
axial form factor Ga- 

There are no terms of lower order. This is in accord with the power counting for the two 
graphs, which predicts a chiral order of for this contribution to Ga- In order to compare 
this with the g^-terms in Ga worked out in [13], one can proceed as follows. Looking at 
fig. ^, and imagining the vector meson lines shrinking to a point vertex (corresponding 
to a limit where the mass My tends to infinity, with GyRy/My fixed), it is intuitively 
clear that the result corresponds to a pion-nucleon loop graph with an 0{(f') contact term 
replacing the vector meson line. Such contributions are parametrized by the two LECs C3 



and C4 in [13], see eq. (5.2). In fact, identifying 

2y/2GyRy 



M2 



-C4, (5.11) 



one reproduces exactly the corresponding terms in the representation based on the pure 
pion-nucleon theory, cf. eq. (^.21) . That eq. ( ^.11 ) is a good guess can be seen like that: 



Comparing the /oA^-coupling used here, namely, the term proportional to Ry in eq. (^.5|), 
to a more conventional one using a vector field representation for the rho field, 

one deduces 

AVlmRy 
QpNN^p = 



-18- 



Using this in eq. ( 5.11 ), we get 

2mMv 4m ' 



C4 



In the last step, we have assumed a universal rho coupling, MyGy = F gpnn 

as well as the KSFR relation My = 2F'^gj^^j^ [20] (see also the recent discussion in the 

framework of effective field theory in ref. [21]). This agrees with the rho-contribution to C4 

found in [16]. Furthermore, there is no rho contribution to the LEG C3 according to this 

work. 

The result of eq. ( 5.11| ) is not surprising for itself, but the agreement of our findings with 



previous resonance saturation analyses demonstrates one very important thing, namely, 
that the variants of the infrared regularization scheme derived in the previous sections are 
consistent with the standard case of infrared regularization in the pion-nucleon sector used 
in [13]. 

As a side remark, we note that the leading order result from the triangle graphs shows 
no t— dependence and therefore gives no contribution to the axial radius. However, the 
leading contribution to the axial radius can also be related to meson resonances, namely, 
to a tree-level exchange of an axial-vector meson. The pertinent calculation can be found 
in [22], where the axial vector meson-couplings to the pions and nucleons are fitted to 
experimental data for GA{t) (for an earlier study based on chiral Lagrangians, see [23]). 
Equivalently, it can be parametrized by a certain LEG, named ^22 in [13] (see eq. (|5.2| )). 
As already mentioned at the beginning of sec. ^, the difference between the two approaches 
just amounts to a resummation of higher order terms. Compared to the leading order 
term, the t-dependent part derived from the triangle graphs is suppressed by factors of the 
small variable a, which is a reflection of the fact that the infrared regularized loop integrals 
preserve the chiral power counting. 



6. Summary 

In this paper we have presented an extension of the infrared regularization scheme that 
allows for an inclusion of explicit (vector and axial-vector) meson resonances in the single- 
nucleon sector of ChPT. For the processes we have considered here, the meson resonances 
do not appear as external particles, and the corresponding power counting rules for the 
internal resonance lines are set up such that the resonance four-momentum is considered 
to be small compared to its mass. The infrared regularization scheme extracts the part 
of the one-loop graphs to which this power counting scheme applies (for any value of the 
dimension parameter d used in dimensional regularization) , while the remaining parts of the 
loop graphs will in general violate the power counting requirements, but can be absorbed 
in a renormalization of the local terms of the effective Lagrangian. 

After a short review of the infrared regularization procedure used for the pion-nucleon 
and the vector meson-pion system in sec. § and ^, respectively, we have combined the 
analyses of these sections in sec. ^. There, we consider the simplest example of a Feynman 
graph where nucleons, pions as well as (vector) meson resonances show up. It is shown how 
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to extract the infrared singular part of such a graph, and the power counting requirements 
are verified. It should be clear from this example how the infrared singular parts of more 
complicated one-loop graphs (with more nucleon, pion and resonance lines) can be worked 
out (some remarks on the general case can be found at the beginning of sec. ^ and in 
sec. 6 of [2]). Finally, in sec. ^, we have applied the extended scheme to compute a vector 
meson induced loop contribution to the axial form factor of the nucleon, and demonstrate 
that the result agrees with the result for Gyi(t) given by Schweizer [13] in combination with 
the resonance saturation analysis for the pion-nucleon LECs in [16]. There are, of course, 
many other possibilities for applications of the scheme developed here. Finally, we add the 
remark that a suggestion for an extension of infrared regularization to the multiloop case 
was made in ref. [24]. 
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A. Explicit expression for 

In eq. (8.4) of ref. [1], a closed expression for ij^y for d — > 4 was given that looks much 
simpler than our result, eq. (|3.7] ), where there is still an infinite sum to be performed. On 
the other hand, it is quite tedious to work out the chiral expansion of the result of [1], due 
to the rather complicated expressions 

2 = ^ + a - 1 ± y^(/3 + a - 1)2 - 4d^^ (A.l) 

used there. These two expressions are nothing but the zeroes of the Feynman parameter 
integral encountered in the computation of the loop integral. The corresponding chiral 
expansions start with 

xi = —a + . . . , 

showing that xi is small and negative, while X2 tends to minus infinity for /? ^ 0+. In 
order to show the equivalence of the two results for I^jy, we employ the following relations, 

(3{xi+X2) = P + a-1, (A.2) 
f3xiX2 = a, (A. 3) 

to write 

(a/3)'= _ 1 (t)'' _ {-irr{2k + m + l) fxi\''^"' 



(1_«_^)2A:+1 (1 + |l)2fc+l ^^^^ m!r(2/c + l) V^2 



-20- 



Inserting this in eq. ( ^ ) yields 

jiR > 2. ^ M^-\a)i-^ ^ ^ T{2k + m + l)T{l - I - fc)r(2fc + 1) / x^\ 
'"''^^^ (4vr)i/3x2 hi^o r(2fc + l)r(fc + l)r(m + l) V xj 



In the second hne we made use of eq. ( A.3 ). Now we change the summation indices 
according to 

j = k + m, I = k = j — m, 



and use the following sum formula for Gamma functions, 

^ r(j + / + i)r(x-0 T{x-m-x) 



for X = 1 — ^. We shall give a short outline of a proof for eq. (A. 5): Dividing this equation 
by r(x — j), both sides are just polynomials in x of degree j, with coefficient 1 in front of 
x^ . Consequently, one only has to prove that both polynomials have the same set of zeroes, 
namely {— 1, — 2, . . . , — j}. This is not difficult, making use of 

^ (-ip-^r(j-^ + i) ^ (-iy-^r(j-z + i) 

^ r(j - 1 + mi + 1) ii^ P> ^ r(j - / + i)r(/ - . + 1) 



E(\')(-ir^-" = (i-ir^ 

n=0 



for < i < J. These remarks should be sufficient to complete the proof of eq. ( |A.5D . 
Returning to eq. (A.4), we employ eq. ([A.5|) to write 



iiR^i^q^) = MtrHa)i~^x, ^ r(j + / + i)r(i - f - /) ^ xi y 



(47r)i j^^U r(j - / + i)r(/ + 1) V 
M^-\a)i-^x, ^ r(f - i)r(i - f - j) /xiy 

(4vr)i r(f-l-j) 1^2^ ■ 



The last line of eq. ( |A.6| ) is exactly the result that was derived in sections 7 and 8 of [1]. 
This can be seen by substituting 

a — > — xi, ^ — X2 , d^- — 2 
in eq. (7.7) of that reference, and multiplying the result with 

r(2 - ^)M^-^ 



-{-PX2)- 

(47r)2 
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as explained at the beginning of sec. 8 of [1] . In the hmit d ^ 4, the series of eq. (^) can 
be summed up to give 



Imv{i\ d^A) = 2xiX - Ixiil - Ino) - (xi - X2) InM - ) ) , (A.7) 



where 



A = Mil ( 1 _ i(in(4vr) - 7 + 1' 
167r2 Vd-4 2^ ^ ^ ' ' 



Eq. (fOP is identical to eq. (8.4) of [1]. 



B. Alternative derivation of I^^^v 

In this appendix, we present an alternative derivation of the infrared singular part of the 
loop integral luBvik'^) (see eqs. ( [4.2| , ^^ )) using the prescription of Ellis and Tang we have 
already explained at the end of sec. ^. However, at some places we will also use Feynman 
parameter integrals, so the derivation outlined here is to some extent a mixture of the 
standard infrared regularization procedure and the method of Ellis and Tang. In complete 
analogy to the steps performed at the end of sec. ^ we start with 

2 f d'^l i {—2k • ly 



{2ttY {{p - ly - m2)(/2 _ M2) ^ (fe2 + /2 _ m2 



dH i ^ {-2k -ly 



{27rY {{p - ly - m2)(/2 _ M2) ^ (A;2 + M2 - M2 

^ f i{-2k • ly 

^ / {2ttY (fc2 + M2 - M2 )i+i((p - If - m^yP - M^)' ^ ' ^ 



3=0 



Now we could use the procedure outlined in [12] to expand the nucleon propagator, together 
with an interchange of summation and integration. In the present case, however, it is easy 
to see that it is equivalent to use the common Feynman parameter trick for the remaining 
loop integrals and extend the parameter integration to infinity like in sec. |2|. For Imb (see 



eqs. (2T) and (|2.3| )), the splitting of eq. (2^) corresponds to 

(B.2) 



{{p-iy -m^){P - M"^) p^ -m? -2l-p + M'^\P -M"^ {p - If - m? 

This has also been noted in [2], see eqs. (22,23) of that reference. On the other hand, the 
first term on the r.h.s of eq. ( p. 21 ) is exactly the integrand that gives the soft momentum 
contribution in the sense of Ellis and Tang, see e.g. eq. (7) in [12]. Therefore, it is consistent 
to continue the series of steps in eq. (|B.1| ) with 

oo 

" ■ ^ ^ ^ (A;2 + M2 - M2 )i+i J {2t:Y [{{p - I? - m'^)z + {P - M2)(l - z)f 
= ir/V(fc'). (B.3) 



dz f dH i{-2k-iy 
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The loop-integration can be done using the following generalization of eq. (|3.8| ): 



n] 



{-ly-^i-k^M^YM"-^'' ^ ; ^ — — f — - (B.4) 



(2vr)'^(/2-M2)'- ^ ' ^ > r(i) (4vr)ir(r) 

for r, n £ N. This gives 



7 



' .t^.^o^2N (7-«-/?y-^ r(i)(4vr)i h i.^-az + a) 



Here we have also used the on-shell kinematics specified in eq. (4.1). The sum over the 



even integers i extends only to j — 1 if j is odd. Defining new indices, 

7 • ^ 7 ^ 

J = 7 , t = — , 

and reordering the series correspondingly gives the following expression for /' 



A/By. 



,=0 (4vr)n7 - « - /3)^+'+ir( J - / + l)r(2/ + l)r(i) Jo (z2 _ a, + 



Using eq. ( p. 10 ), it is straightforward to see that this equals I^jby of ^1- (Eil)) expected 



(the remaining parameter integral can be done with the help of eq. ( |4.8| )). 
C. Loop Integrals 

Here we list the decomposition of the loop integrals with tensor structures in the numerator, 
which we need in sec. ^. All loop integrals in this appendix are understood as the infrared 
singular parts of the full loop integrals, but we will suppress the superscript IR for brevity. 
As a consequence, all loop integrals that do not contain a pion propagator are already 
dropped here, since they have no infrared singular part. Also, we will use the mass shell 
condition = for the nucleon momentum p. We start with 



where 



In complete analogy. 



(27r)'^((p-/)2-m2)(/2-M2) ^ 



p'I'm'b^ (C.l) 



/ 



(27r)'^((/-A:)2-M2)(;2_M2) 
with 
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Integrals of type MB and MV are also needed with a tensor structure l^l'^ in the numerator. 
They are decomposed as 



(27r)'i((Z-/t)2-M2)(/2_M2) 



(C.3) 



where the coefficients of the tensor structures are given by 



{d-l)tflyik) 
{d-l)t%{k) 



-(AfV H TTo -'M, 



4A;2 

ttt; Imv 7T-^ 



4A;2 

U2 I /\#2 



4fc2 



4fc2 



The corresponding coefficients in the meson-baryon case, t^MB(p)^ derived from 

these results by substituting k ^ p, My m. 

We turn now to loop integrals with three propagators. First the vector integral: 



{2-kY Hp - I? - m^){{l - - M2 )(/2 _ M2) 



'+PYlUv + {k-pYlfjj,y, (C.4) 



with 



1 



^MBV 2fc2(4m2 - A:2) 

- k^lMV - {2rr? - k^)lMB 



{2m' + (2m2 - k'^k' - M{.))Imbv 



^MBV 



2fc2(4m2 - A:2) 
+ 2,k'^lMV - {k^ + 2m^)lMB 



{2M\m' - k') + {k' - M{,){k' + 2m'))l 



MBV 



We remind the reader that we use = m? = {p — k)"^ here. The scalar loop integral with 
three propagators occuring in this decomposition was named Imbv in sec. eq. (^). 
However, we noted there that it is equal to Imbv for on-shell nucleon momenta. 
The tensor integral 



^MBV 



dH 



{27tY {{p - /)2 - m2)((/ - A;)2 - M2 )(/2 _ M^) 



(C.5) 



can be decomposed as 



i^Bv = 9^"Ci+ip+kriP+krc2Hp-kr{p-krc3Hip+krik-prHp+krik-pr^^ 
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Here the coefficients Cj are given by 

1 - k'^ 
M^Imbv --{k^-M^ + 2M^)lUv + ^ Imbv 



m 



'^M^Imbv + 



^2(4^2 _ ^2) ^ 

4 ^MB - X^mV -m {d- l)Ci 

C3 = p(4^2_fc2) ^'("*' + 2fe')/MBy - ^(fc' -M^ + 2M2)((m2 + 2A;2)/j^5^ 

2vr? + ^(1) 



P(4m2 — P) 



+ (^/e - m )1mbv j + ^-'My 4 
- {m'^ - k'^){d - l)Ci 
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